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Abstract 

The discrete energy-eigenvalues of two nucleons interacting with a finite-range nuclear force and 
confined to a harmonic potential are used to numerically reconstruct the free-space scattering 
phase shifts. The extracted phase shifts are compared to those obtained from the exact continuum 
scattering solution and agree within the uncertainties of the calculations. Our results suggest that 
it might be possible to determine the amplitudes for the scattering of complex systems, such as 
nd, nt or na, from the energy-eigenvalues confined to finite volumes using ab-initio bound-state 
techniques. 
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I. INTRODUCTION 



Quantum scattering of strongly interacting few-nucleon systems is complicated and requires 
careful treatment of asymptotics, antisymmetrization effects as well as the dynamics gener- 
ated by nuclear forces. Full treatments of antisymmetrization with correlations have become 
routine in bound-state and quasi-bound-state solutions of light nuclei using ab-initio tech- 
niques based on nucleon-nucleon (NN) and three-nucleon interactions (3A^). Techniques, 
such as the No Core Shell Model (NCSM) (see, e.g., Ref. P E]), Green's Function Monte 
Carlo (GFMC) (see, e.g., Ref. [3l H]), and the Coupled Cluster (CC) approach (see, e.g., 
Ref. [3 E]), are used to calculate ground and excited states of light nuclei. The precision of 
these calculations has reached a point where further progress is now limited by the fidelity 
of the input interactions. 

In reaction calculations of scattering properties of light nuclei, progress has been less 
pronounced, though nonetheless significant. R-matrix analysis (see, e.g., Ref. [TJ |S]) has 
been historically the empirical workhorse, providing impressive fits to a range of experimen- 
tal data. More microscopic approaches to the scattering of light nuclei are based on the 
Resonating Group Method (RGM) |9]-[TT]. A promising avenue for ab-initio calculations of 
scattering of light ions comes from the coupling of the RGM reaction method with the NCSM 
bound-state technique [12] . For this method, the large computational resources required to 
achieve convergence provide the limiting constraint on reliably calculating scattering param- 
eters for processes with A > 5. It would be significant if existing bound-state techniques, 
and their accompanying precision, could be further exploited to reliably determine scattering 
amplitudes for multi-nucleon systems. 

During the last twenty years the general technique of effective field theory (EFT) has 
been developed and applied to multi-nucleon systems. Effective Field Theory provides a 
description of observables, consistent with the approximate chiral symmetries of quantum 
chromo dynamics (QCD), in terms of a small number of expansion parameters within a 
plane-wave basis. These expansion parameters are used to relate the experimentally de- 
termined scattering parameters and bound-state properties of few-nucleon systems to the 
coefficients of operators at a given order in the EFT expansion. In principle, this allows 
for systematically improvable calculations of multi-nucleon observables. Applying the EFT 
framework within an oscillator basis has also been investigated [T314T5] . More recently, it 
has been suggested that the EFT framework might be fruitfully applied to multi-fermion 
systems confined in a harmonic potential [T6H18], and might be usefully married with the 
NCSM calculational scheme. In Ref. [19] it was suggested that the scattering properties of 
certain complex nuclear systems could be calculated from the spectrum of the same systems 
confined to a harmonic potential. This was demonstrated for two confined particles at the 
unitary limit in Ref. [20] • A lattice formulation of EFT coupled with an external harmonic 
potential is currently being developed [21j . 

In this work we investigate the simple two-nucleon system confined in a harmonic poten- 
tial of the form Vho = ^M^u^'^r'^, and interacting via nuclear forces in uncoupled partial 
waves. Since two-body techniques are well established for both scattering and bound states, 
this system is ideal for determining the extent to which continuum scattering amplitudes can 
be recovered from bound-state information. An analytic expression that relates the eigen- 
values of two interacting particles moving in a harmonic potential to the scattering phase 
shift at those energies, analogous to "Liischers method" [22142^ that is used in Lattice 
QCD, allows for the scattering phase shift to be determined in the limit that the oscillator 
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length is large compared to the range of nuclear forces. As the latter is characterized by the 
Compton wavelength of the pion, with increasing confinement this leads to modifications to 
nuclear forces due to the harmonic potential that must be systematically removed in order 
to accurately predict the scattering amplitude. This is achieved by calculating the energy- 
eigenvalues over a range of harmonic frequencies, u, determining the scattering amplitude 
over this range of u, and then extrapolating to the u = limit. We verify the extracted 
phase shifts by comparing them to the results of an independent scattering calculation. 
The two methods are found to yield the same phase shifts within the uncertainties of the 
calculations. 

II. PHASE SHIFTS FROM THE ZERO-RANGE RELATION 

It is well-established that the energy-eigenvalues of an interacting system of particles confined 
to a finite volume can be used to determine the scattering phase shift (at the energy- 
eigenvalues) when the size of the confining region is much larger than the range of the 
interactions between the particles. For instance, relating the scattering phase shift to the 
energy-eigenvalues of two-nucleons confined to a spherical region by solving the Schrodinger 
equation with a Dirichlet boundary condition is a problem that appears in standard texts 
on nuclear physics, see, e.g., Ref. |25jj. This method has been successfully employed in the 
latticization of low-energy EFT's to predict (to a given level of precision) the ground state 
energies of light nuclei and their volume dependence [26] . 

Volume dependence is also used to determine meson-meson, meson-baryon and baryon- 
baryon elastic scattering phase shifts from the energy-eigenvalues of these systems calcu- 
lated with Lattice QCD ^ (for a recent review, see Ref. [2H])- Lattice QCD calculations are 
generally performed in spatial volumes with cubic symmetry and with periodic boundary- 
conditions (BC's) imposed upon the fields at the edges. This reduces the number of momen- 
tum modes, and hence reduces the kinetic-energy contributions to the calculated processes. 
The non-relativistic relation between the energy-eigenvalues and the scattering amplitude 
(the "Liischer relation" ) has been shown to be valid even in quantum field theory [231 El] ■ 
Given the energy-splitting, AEn, between the two-hadron state and the hadron masses, m, 
the real part of the inverse scattering amplitude below inelastic thresholds is 

P„cot5(p„) = ^S^^j j , S(^) ^ - , (1) 

where AEn = '^\/ v\ + ~ 2m, 5(pn) is the energy-dependent phase-shift, and the 
limit A — 7- oo is implicit. The S-function is the Green function, G//iy(0,0), for two-hadron 
plane- wave eigenstates (and straightforwardly generalizes to hadrons with different masses). 
The Liischer relation between the scattering amplitude and the energy-eigenvalues in the 
finite lattice volume, given by Eq. ([T]), is valid when the spatial extent of the lattice is large 
compared to the range of the interaction, R. Corrections to the relation are found to behave 
as ~ e~^/'^, see, e.g., Ref. [22] ■ If, in the continuum (infinite- volume) limit, a system contains 
a shallow bound state, as is the case in the NN ^S\-^Dx coupled channel, the periodic BC 



The Maiani- Testa theorem [37] precludes the extraction of scattering matrix elements from Euclidean- 
space Green functions in the infinite- volume limit except at kinematic thresholds. 
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in a finite volume increases the binding energy of tlie state. Tlie finite-volume corrections 
scale as e~"'°^ [30], where 70 is the binding momentum in the continuum. In contrast, the 
continuum scattering states have power-law dependences upon the lattice extent for L 3> -R, 
with energies that behave as ~ for the ground state and ~ for the higher-energy 
states. 

The systems that we considered in this work are comprised of two nucleons in a harmonic 
potential (with oscillator frequency u) interacting through NN forces. We will consider 
the JISP16 potential |3TH33] . which reproduces the low-energy NN scattering data with 
a x^/dof ~ 1.0, but our results are general and the technique can be applied to other 
NN interactions. The EFT-method that was used to (re-)derive the Liischer relation in 
Eq. ([1]) in non-relativistic quantum mechanics |30], can be used to (re-)derive the relation 
between p^^^^ cot 5 i{p) in an uncoupled partial wave with angular momentum I and the 
energy-eigenvalues of two nucleons in a harmonic potential , 

p^'+^cot^Kp) = (-1)'+^ {2mur"' ^ ^ _ ! , (2) 

where e = E/u and E = /m is the fully interacting energy in the center of mass frame. 
While the EFT derivation using the pionless-EFT is valid only up to momenta associated 
with the cut in the t-channel from the exchange of one pion, the relation is valid up to the 
inelastic threshold. Equation ([2]), like the Liischer formula, is valid in the limit of zero-range 
interactions. The harmonic potential, being non-zero everywhere, except at the origin, 
modifies the interaction between the two nucleons, and the A^A^ phase shift at the outer 
range of the nuclear potential differs from that in free space. This is a finite-range effect, 
and unlike the situation encountered in Lattice QCD calculations, it is not expected to be 
exponentially suppressed (in co). Equation ([2]), in conjunction with the leading order (LO) 
term in the effective range expansion (ERE) of p cot 6, the scattering length a, has been 
used to determine the spectrum of dilute cold atoms in traps with essentially zero-range 
interactions, particularly in the vicinity of Feshbach resonances [351438] . On the other hand, 
smce Eq. relates the energy-dependent phase shift to the energy-eigenvalues E of the 
confined system, knowledge of the spectrum of the two-particle system allows the extraction 
of the continuum scattering amplitude up to finite-range corrections. 



A. Loosely Bound States in Weak Harmonic Potentials 

For attractive S'-wave interactions with positive scattering length, > 0, a bound state 
exists and the binding energy Bq can be written in terms of the binding momentum, 70, 

m 

where 7 ~ 1/a for scattering lengths large compared to the range of the interaction. Refining 
the estimate of the binding energy gives 7 as the solution to 

^ + Jro7o - 70 = , (3) 

where tq is the effective range of the interaction, and the ERE of pcot 6 = — ^ -|- ^rop^ -|- ... 
has been truncated at second order, which, in the case of 5'-wave interactions between 
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nucleons, is sufficient for most purposes. The presence of the harmonic potential, and in 
particular its non-zero value throughout the volume of the bound-state, gives rise to a power- 
law modification to the binding energy [33] , even in the limit of zero- range interactions. The 
location of the state corresponding to the free-space bound-state can be found directly from 
Eq. g in the zero-range limit, and for small u the shift in the energy of the bound-state is 
perturbative in u"^, 

= Bo - r ^ + Oiu^) = Bo - CzRUj' + Oiu^) , (4) 

8(1- 7oro) Bo 



where we define Czr = [8Bo (1 — 7o?"o)]~ for later reference. Equation (|4|) indicates that, 
given the bound-state energy B^^ calculated at different values of u, the continuum binding 
energy Bo could be determined by an extrapolation in to = 0. This same extrapolation 
can also be done in the presence of finite-range corrections since, as we show later, these 
corrections occur at order u"^ for small u. In the ^Si-^Di coupled-channels that contain the 
deuteron with binding-energy Bo = 2.224575 MeV (70 ~ 45.7 MeV) and with an S-wave 
effective range of tq ~ 1.74 fm, the coefficient is (7^^^=0.0944 MeV~^. 

The LO shift in the bound-state energy given in eq. Q can be recovered from the bound- 
state wavefunction based on ERE, 

^(^^Hr) = ^sho.t(r) + .[f ^ — , (5) 

V 27r VI - loTo r 

where ipshonif) is the short-distance component of the wavefunction that has support over 
a radius r <ti . The factor of 1/a/1 — 7oro in eq. (pi) is determined by the residue of 



the pole in the scattering amplitude. At LO in perturbation theory, the contribution to the 
energy of this state from the harmonic potential is 

1 1 o;^ 

AEo = (V^^^^^I-mwVV^'^^^) = —, ^ IT- + short - distance , (6) 

2 8(1- 7oro) Bo 

in agreement with the result in eq. (|4]). 



B. Scattering States in Weak Harmonic Potentials 

It is useful to construct perturbative expansions for the energy-eigenvalues in the zero-range 
limit. As we show later in sect. IV A these relations can be used to readily extract effective 
range parameters given the low energy spectrum of the system. Using the zero-range relation 
given in eq. ([2]) it is straightforward to determine the location of the energy-eigenstates in 
the limit that y/muj/{p cot 5) ^ 1, and also in the unitary-limit where y/muj/{p cot 5) ^ 1. 
In the y/muj/ {p cot 6) ^ 1 limit, the g*'' energy-level with orbital angular momentum / is 
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where Eq = (| + Z + 2g) with b = l/y/mu, and the H{x) are harmonic numbers 
When apphed to the lowest-lying S-wave state, one finds that, for small-w. 
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where p cot(5o is evaluated at p'^/m = 3u/2 ^. The first excited S-wave state is located at 
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3(6 log 2 - 5) 3 (37r2 - 4(11 + 9 log 2(3 log 2 - 5))) 



bp cot^o vv^27r 47r bp cotSo 



6V27r3/2 (bp cot5o) 



+ . . . 
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where p cot(5o is evaluated at p^/m = 7u)/2. The finite range corrections to these expressions 
can be introduced by replacing pcot^o pcot^o + ^^'^ + ••• for small u, where the finite- 
range corrections depend upon the interaction, and cannot be determined from scattering 
parameters alone. It makes little sense to continue the expansion in y/moj / [p cot 6) to higher 
orders due to the appearance of the range corrections. 

The same expansion can be applied to the P-waves, for which the expansion is in terms 
of (ma;)^/^/(p^cot5i) ^ 1. The lowest-lying continuum state is located at 



E, 



(1) 



6(31og2-4) 
2 ~ 63p3cot(5i ~ TT (6Vcot5i)^ 



+ . . . 



(10) 



with p^cot^i evaluated at p'^/m = 5ijj/2, the first excited state is located at 
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(1) 



2^ 



15(30 log 2 -31) 



2 6Vcot(5i An (6Vcot5i) 



+ . . . 



(11) 



^ Our definition of the oscillator parameter differs from that of Ref. (20j by a factor of y/2. 
To recover the results of Ref. [20], p cot(5o is replaced by the ERE evaluated at p'^/m — 3w/2, 

1 3 To 

p cotdo = h-— ma; + . . . , 

flo z z 

and eq. (|8|) is then re- arranged in powers of b^^. 
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with p^cot^i evaluated at p^/m = 9a;/2, and the second excited state is located at 



E'i'^ 13 105(4101og2 - 389) 



UJ 2 6Vcot(5i 12877 (feVcot^i)^ + • • • , (12) 

with p^cot(5i evaluated at /m = 13u/2. 

The limit of large scattering length, a/b ^ 1, and small range, r/b <^ 1, the unitary limit, 
can also be considered. An expansion in powers of pcot5Q/{muy^^ can be performed, and 
the lowest-lying S-wave state is located (by expanding about the poles in the denominator 
of eq. ) at 

1 /2"pcot(5o 



where pcot^o is evaluated at p'^/m = u/2. This generalizes the results of Ref. [35] . 



III. A TOY MODEL 

For a harmonic potential with an arbitrary value of u one must rely on eq. ^ to extract con- 
tinuum phase shifts and scattering parameters from the location of the energy-eigenvalues, 
while keeping in mind that finite-range effects are present and will move the calculated phase 
shift away from its true value. To test the utility of eq. ^ in the presence of a finite-range 
interaction and to develop a "feel" for the size of the finite-range corrections, we use the toy 
example of two particles interacting via a spherical well. To make this system as 'nuclear- 
like' as possible, the depth and width of the well are tuned to reproduce gross features of 
the deuteron system. In particular, with a well depth Vq = 48 MeV and radius Ro=1.7 fm, 
the single bound state has a binding energy i?o=2.22 MeV. The scattering phase shift for 
this potential is known to be 



6o = tan 



-1 



tan (^/Ro^I (^lab + 2Vo] 



E,.^RlVofi , (14) 



where fi is the reduced mass. 

This system is placed within a harmonic potential and the resulting two-body spectrum 
at various oscillator frequencies is determined numerically. For each oscillator frequency, 
the spectrum is used to extract the scattering phase shift by virtue of eq. ([2]). Because 
the spectrum is discretized, the extracted phase shifts occur at discrete points. By varying 
the oscillator frequency, the energies at which the phase shift is determined vary thereby 
allowing for the energy-dependence of the phase shift to be mapped out. 

For modest-sized oscillator frequencies {uj < 4 MeV) the extracted phase shifts agree well 
with the exact result given in eq. (14) (within 0.1%), as shown in Fig. [T| as the effects of 



the harmonic potential are negligible within the range of the spherical well. The situation 
changes, however, for large oscillator frequencies, also shown in Fig. [T] In this case the exact 
phase shifts and extracted phase shifts have appreciable differences due to the finite range 
of the spherical well. Not surprisingly, the confining nature of these potentials distorts the 
interaction of the two particles within the spherical well, which is demonstrated in Fig. |2} 
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FIG. 1: (Color online) Extracted phase-shifts for the spherical- well toy-model for oscillator fre- 
quencies from uj = 0.5 MeV to w = 4.0 MeV (left panel), and from uj = 8.0 MeV to w = 64.0 MeV 
(right panel). For each oscillator frequency, the phase shift was determined from the lowest eleven 



energy-eigenvalues (excluding the bound state). The exact continuum phase shift, given by eq. (14), 
is the solid black curve. Appreciable deviations in the phase shift at larger oscillator frequencies 
are due to the finite range of the spherical well. 




r(fm) 

FIG. 2: (Color online) The potential between two particles interacting via a spherical well, Vgw^ 
confined by a harmonic potential, Vho^ for different oscillator frequencies. As the oscillator fre- 
quency increases the distortion of the spherical well increases. 

An interesting feature of these finite-range effects is that for a given oscillator frequency, the 
effects are largest at lower energy, and diminish as the energy of the system increases ^. 



^ In the high-energy limit, in which the nucleon wavelength inside the range of the nuclear interaction is 
small compared to the length scale over which the potential varies significantly, the LO contribution of 
the harmonic potential to the s-wave phase shift, calculated in the WKB approximation, is 



2V2 



dx X 



^E-Vnn{x) 



E 



(15) 



where is the reduced mass of the two- nucleon system, and Vnn{x) is the (central) NN potential. In 
the case of a toy model of NN interactions where a spherical well of depth Vq and radius Rq is used to 
describe the NN potential (Vo and _Ro are tuned to reproduce the scattering length and effective range). 
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IV. REALISTIC NUCLEAR FORCES 



It is important to determine how well this method works for realistic NN interactions. There 
are a number of modern NN potentials that could be used for this numerical comparison, but 
for simplicity we work with the JISP16 potential pTH55] . This A^A^ interaction is constructed 
so as to reproduce the measured NN scattering phase shifts to high precision over a wide 
range of energies, below -Eiab^ 350 MeV, and is known to provide a good description of p 
shell nuclei |311 [33] without an additional 3A^ interaction. It was developed using inverse 
scattering techniques, followed by off-shell tuning with phase-equivalent transformations to 
describe selected light nuclear properties up to ^^O. Using this interaction, the spectrum of 
the two confined particles was found by diagonalizing the Hamiltonian in the relative HO 
basis space for each partial wave. The size of the HO basis was increased until the spectrum 
converged to a prescribed precision. In order to access the lower energy phase shifts, we 
decreased u which consequently required an increase in the size of the basis space to obtain 
convergence. This limited the range of small u that we investigated (w > 0.4 MeV with a 
maximum basis dimension of 1800x1800). In order to improve convergence with increasing 
basis-space dimension, the choice of the HO frequency for the basis space was adjusted 
independently of the frequency of the external confining potential. 

In Figs. |3]to|6]we show the application of eq. ([2]) to four different partial waves in the A^A^ 
system: (/ = 0), ^Pq {1 = 1), ^£'2 {I = 2) and ^Fg (/ = 3). The extracted phase-shifts 
were obtained from the low-lying spectrum of the A^A^ system in harmonic potentials with a 
range of frequencies (the points in each figure). For comparison, the phase-shifts calculated 
by solving the Schrodinger equation in the absence of the harmonic potential are shown as 
the solid curves in each figure. 




50 100 150 200 250 300 350 5 10 15 20 25 30 



£iab(MeV) £iab(MeV) 

FIG. 3: (Color online) The phase-shift in the ^5o-channel, 5aj(-E^iab)j evaluated at the energy- 
eigenvalues determined for a range of values of uj defining the external harmonic potential, using 
eq. The solid curve corresponds to the phase-shift, 5a;=o(£'iab)) determined by a direct evalua- 
tion in free-space. The right panel is a magnification of the left panel. 



the correction to the phase shift is 

1 / /; \ 3/2 r°° 1 / X 3/2 

SUE) ~ S^ME) ^ {'^) io'j^ dx V^j,{x) = ^ ['^) u;'VoRl ■ (16) 
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FIG. 4: (Color online) The phase-shift in the ^Po-channel, (5(^(-E'iab)! evaluated at the energy- 
eigenvalues determined for a range of values of oj defining the external harmonic potential, using 
eq. The sohd curve corresponds to the phase-shift, 5^=Q{E\.^\y), determined by a direct evalua- 
tion in free-space. The right panel is a magnification of the left panel. 
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FIG. 5: (Color online) The phase-shift in the ^D2-channel, (^tj(-E'iab)) evaluated at the energy- 
eigenvalues determined for a range of values of w defining the external harmonic potential, using 
eq. (U. The sohd curve corresponds to the phase-shift, (^aj=o(-E'iab)) determined by a direct evalua- 
tion in free-space. The right panel is a magnification of the left panel. 



A. Numerical Analysis 

The harmonic potential modifies the interactions between the two nucleons at all distance 
scales, and as such, there are modifications to the potential between the nucleons over the 
range of the nuclear forces, leading to short- distance corrections to the relation between 
p cot b and the energy-eigenvalues given by eq. (|2]) . The energy-eigenvalues are calculated 
in a given energy-interval for a range of values of lo in order to extrapolate the phase- 
shift 5aj(-Eiab), to the a; = limit, (5a;=o(-^iab)5 and hence eliminate the modifications to the 
nuclear force due to the harmonic potential. This procedure is not as straightforward as it 
naively appears due to the fact that for each value of the energy-eigenvalues (generally) 
have different values, and an interpolation of (5(^(-E'iab) within the energy-interval is required 
for each uj in order to extrapolate to 5tj=o(-E'iab) at any given energy ^ . For the sake of 



^ For a given energy, a range of values of uj could be iteratively tuned to produce the same energy-eigenvalue. 
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FIG. 6: (Color online) The phase-shift in the ^i<3-channel, 6ui{Eii^\,), evaluated at the energy- 
eigenvalues determined for a range of values of uj defining the external harmonic potential, using 
eq. The sohd curve corresponds to the phase-shift, (5a;=o(-E'iab); determined by a direct eval- 
uation in free-space. In the left panel this solid curve is obscured by the points calculated using 
eq. Q. The right panel is a magnification of the left panel. 

demonstration, we focus on the phase-shift in the ^So channel, but the methodology can be 
applied in all channels. 

The energy-eigenvalues of two nucleons interacting in the ^5*0 in a harmonic potential 
were calculated for a range of values of u from u = 0.4 MeV to u = 15.0 MeV. For each 
eigenvalue, the scattering phase-shift 5tj(-E'iab) was calculated using the zero-range relation 
in eq. ([2]), the results of which are shown in Fig. |3} The "exact" phase-shift Sco=o{Eiah) that 
is determined by solving the Schrodinger equation for the phase-shift in the absence of the 
external harmonic potential is shown in Fig. [s] as the solid curve. For ^ 1.0 MeV, and 
for the energy-eigenvalues shown in Table [1} the phase-shift calculated from the zero-range 
relation in eq. ^ is very close to the actual phase-shift. For larger values of u there are 
noticeable deviations from the exact result, but these deviations are found to become smaller 
as the energy increases. 

The energy- dependent interpolation of the phase-shift for a given u that is required in 
order to extrapolate (5cj(-E'iab) to (5tj=o(-E'iab) is the most problematic part of this numerical 
analysis. In the range of energies for which the ERE is formally convergent (|p| < m7r/2 
resulting from the location of the t-channel cut in the one-pion exchange amplitude) an 
expansion of p cot 6 in powers of the energy reproduces the scattering amplitude. In contrast, 
for energies outside of this range, but below the inelastic threshold (for which the relation 
between p cot S and the energy-eigenvalues in the harmonic potential in eq. ^ remains valid) 
such a power-series does not describe the amplitude. As such, without directly solving for 
the amplitude, one does not know the form to use for the interpolation in relative energy 
E beyond |p| = We do not attempt to resolve this issue, and restrict ourselves to 

the energy range for which the ERE formally converges ^. In Fig. [7]we show the extracted 
values of p cot 5 as a function of relative energy E for a; < 1 MeV. For each a; < 1 MeV a 
fourth-order polynomial in E is fit to the values of p cot S shown in Fig. [T] (the order was 
chosen to minimize the x^/dof of the fit and to achieve a stable fit under the change of 



^ Within this range, this part of our analysis is formally equivalent to the pionless EFT description given 
in Refs. [UlEl]. 
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TABLE I: The lowest eight energy-eigenvalues in the center-of-mass frame and their associated 



phase-shifts found from eq.(|2|) in the So-channel for uj < 1 MeV. 





CO = 0.4 MeV 


oj = 0.5 MeV 


u = 0.6 MeV 


CO = 0.8 MeV 


UJ = 0.9 MeV 


u = 1.0 MeV 


El 


0.66642 


0.81618 


0.96488 


1.2610 


1.40898 


1.55711 




58.5279 


60.0449 


61.1382 


62.5816 


63.0684 


63.4511 


E2 


2.22732 


2.78198 


3.33893 


4.4597 


5.02344 


5.58933 


6UE2) 


64.2586 


64.45411 


64.4124 


63.9721 


63.6576 


63.3088 


E3 


3.82836 


4.7907 


5.75672 


7.69914 


8.67518 


9.65424 


Su^iEs) 


64.2495 


63.7768 


63.1847 


61.8606 


61.1775 


60.4948 


Ea 


5.4363 


6.80548 


8.17932 


10.9398 


12.326 


13.7158 


5UE4) 


63.3856 


62.4755 


61.5194 


59.6034 


58.6685 


57.754 


E5 


7.04598 


8.82126 


10.602 


14.1783 


15.9733 


17.7726 


SUE5) 


62.3082 


61.0673 


59.832 


57.4512 


56.3152 


55.214 


Eq 


8.65604 


10.8369 


13.0238 


17.4144 


19.6173 


21.825 


SUEe) 


61.1822 


59.6694 


58.2013 


55.4272 


54.1198 


52.8615 


Ef 


10.266 


12.852 


15.4446 


20.6482 


23.2584 


25.8736 


SUEr) 


60.0611 


58.3139 


56.6418 


53.5231 


52.0696 


50.6769 


Es 


11.8758 


14.8665 


17.8645 


23.8802 


26.897 


29.9188 


SUEs) 


58.9638 


57.0076 


55.1527 


51.7288 


50.1421 


48.6454 
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FIG. 7: (Color online) pcot6 as a function of the energy in the center-of-mass frame in the ^Sq- 
channel extracted from the energy-eigenvalues of the two-nucleon system in harmonic potentials 
for a range of oscillator frequencies lo. 

order ^ ). With the interpolating functions, it is then possible to choose a particular value 
of E and extrapolate p cot S^{E) to p cot 6i^=o{E) , from which the phase-shift 6^=o{E) can 
be recovered. The ^-extrapolations at E = 1 MeV and E = 5 MeV are shown in Fig. |8} A 
fit function of the form pcot 6^^ = A + B u"^ is used to extrapolate to cu = 0, as also shown in 
Fig. |8| The small observed scatter of the points about the best fit function is attributed to 
the form of the interpolation in E (and the increasing separation between energy-eigenvalues 



^ A full systematic study of uncertainties would include the variation of the fit with the polynomial order. 



12 



0.07560 p 
0.07555 - 



,S 0.07550 - 



3 0.07545 - 
O 0.07540 - 

o 

0.07535 - 




(jj^ (MeV^) 



0.19850 - 



3 
=0 



Oh 



0.19845 - 



0.19840 L 




0.0 0.2 0.4 0.6 0.8 



co^ (MeV^) 



FIG. 8: (Color online) The extrapolation oi pcoi5^{E) to pcoi5u}=o{E) in the ^S'o-channel at 
E = 1 MeV (left panel) and E = 5 MeV (right panel). The solid lines correspond to the best fits 
of the form pcot 6^} = A + B u'^, and the dashed lines correspond to the 99% confidence intervals. 
The red points with their associated 1-sigma uncertainties correspond to 5uj=o{E). 



with increasing u), and not the finite model-space as the energy-eigenvalues have converged 
to six significant digits which we establish by increasing Nms sufficiently. An important point 
to note is that the results of the calculations at the smallest few values of u are all within 
~ 0.1% of the extrapolated values. Therefore, to determine the phase-shift at this level of 
precision, no extrapolation in u"^ is required. The extrapolated phase-shift 5a;=o(-£'iab) in the 
^S'o-channel is shown in Fig. [oj and is found to agree with the exact phase-shift (the solid 
curve) within the uncertainties of the calculation. The points with uncertainties correspond 
to the phase-shift derived from the energy-eigenvalues extrapolated to w = evaluated 
at regular intervals in E. Uncertainties in the extrapolated phase-shifts, which are at the 
~ 10~^-level, can, in principle, be reduced further by calculating at even smaller values of 
u. 

We have numerically explored some of the higher partial-waves. The methodology in the 
higher partial- waves is the same as in the ^S'o-channel. The harmonic potential modifications 
to the nuclear force are seen to increase with increasing partial-wave. This behavior is 
expected due to the fact that the centripetal barrier, and the associated behavior of the 
wavefunction near the origin, forces the wavefunction to larger values of r (but within the 
range of the nuclear force) and hence to larger values of the harmonic potential. Calculations 
at smaller values of u than employed for the S-wave case are required in order to achieve the 
same level of precision, consistent with the conclusions of Ref. [37j. The extracted values of 



the phase-shift in the Pq, D2 and F3 channels extrapolated to a; = are shown in Fig. 10 



In all channels, the extrapolated phase-shifts are found to agree with the "exact" phase-shift 
within the uncertainties of the calculation. 

Determining the energy levels of two nucleons in a harmonic potential involves calculating 
the matrix elements of the full Hamiltonian, including the harmonic potential, in a large 
model-space, with a cutoff on relative excitation energies denoted by ooNus- In the limit 
that Nms 00 the energy-eigenvalues found by diagonalizing the Nms x Nms Hamiltonian 
will coincide with the exact energy-eigenvalues. For a finite-dimensional space, the energy- 



eigenvalues deviate from their infinite model-space values as shown, for instance, in Fig. 11 



making the quantification of the convergence of eigenvalues with respect to A^ms highly non- 
trivial. We do not attempt to resolve this issue here, and all of the energy-eigenvalues we 
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FIG. 9: (Color online) The w-extrapolated phase-shift (5(^=o(-E'iab) as a function of ii^iab in the 
^5o-channel. The points and their uncertainty are determined, at uniform intervals in i^iab, by the 
interpolations and extrapolations described in the text. The solid curve corresponds to the "exact" 
phase-shift. The inset is a magnification around -Biab = 7 MeV that shows the precision of the 
calculation. 




FIG. 10: (Color online) The w-extrapolated phase-shift (5a;=o(-E'iab) as a function of E'lab in the 
^Po-channel (left panel), the ^D2-channel (center panel) and the ^Fa-channel (right panel). The 
insets are a magnification around -Eiab = 7 MeV that shows the precision of the calculations. 

have used in this work have converged to at least six significant digits. 

In this work we have only analyzed uncoupled channels for simplicity. In general, due 
to the spin of the nucleon, and the fact that two nucleons can have 5 = 1, many two- 
nucleon states with total angular momentum J are a linear combination of two orbital 
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FIG. 11: (Color online) The 20**^ energy-eigenvalue in the center-of-mass frame for the ^5o-channel 
with uj = 0.4 MeV as a function of the inverse cut-off of the model-space, 1/Nms- 



angular momentum states, such as the ^Si-^Di coupled-channel which contains the deuteron. 
The zero-range relation between the energy-eigenvalues in the harmonic potential and the 
scattering phase shift, given in eq. ([2]), will be modified to be a relation involving the three 
scattering parameters that describe a two-component coupled-channels system, e.g. 6o, 62 
and 61, and not just a simple relation between one phase shift and the energy-eigenvalues. 
Such relations remain to be constructed for two nucleons in a harmonic potential. As the 
^Si-^Di coupled-channels system contains the deuteron as a bound-state, we can explore the 
behavior of the lowest energy-eigenvalue as a function of u. The perturbative corrections to 
the location of such bound states due to the presence of the harmonic potential are given 
in eq. Q, where the LO deviations scale as ~ 00'^. The binding energies are found to be 
£^0 = -2.2209 MeV, -2.2163 MeV, -2.2098 MeV, and -2.2017 MeV in harmonic potentials 
with oj = 0.2 MeV, 0.3 MeV, 0.4 MeV and 0.5 MeV, respectively. The extrapolation of these 
values to u 



En 



is shown in Fig. 12 
+ D + F CO 



The results are fit well by a polynomial of the form 



where C,D,F and E, 



(a;=0) 



are fit variables, for the 



range of co for which the calculations have been performed. The deuteron binding energy 
extracted from the extrapolation to a; = is eI^~'^^ = —2.22466(4) MeV (which is to be 
compared with the input value of —Bq = —2.224575 MeV). The coefficient of the o;^ term 
is Cfit = 0.0939(4), which is consistent with the value expected in the zero-range limit of 
CzR = 0.0944 from eq. (|6]). One expects both the LO short-range contributions from w 7^ 
and the small D-state admixture due to the tensor force to also depend upon co^, and to 
modify the value of C away from Czr, but it is clear from this work that such deviations 
are small. 

By looking at different energy eigenvalues, the effective range parameters can be extracted 
through the relation 



p2'+^cot5i(p) = -1/a/ + 1/2 + 



-1/ai + 1/2 rimE + 



(17) 



where E is any relative energy eigenvalue which is low enough to ensure convergence of 
the ERE. For example, the low-energy spectrum of the confined ^5*0 system can be used 
to extract the scattering length and effective range using eq. ([2]). We show these extracted 



parameters at the 1-a level in Fig. (13). These extracted parameters vary with u in a 



way that is consistent with expectations and converge to the exact result. For a system 



with ERE parameters that are of natural size, the perturbative expressions from sect. II B 
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FIG. 12: (Color online) The deuteron binding energy as a function of uj"^. The solid line corresponds 
to the best fit of the form Eq = E^~'^^ + C co'^ + D uj"^ + F u:^ , and the dashed lines (practically 
indistinguishable from the solid line) denote the 68% confidence interval. The red point corresponds 
to the ground state energy obtained by extrapolating to a; = 0. The uncertainty is within the size 
of the red point. 
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FIG. 13: (Color online) The extracted ^Sq (left panel) and ^Pq (right panel) scattering lengths 
and effective ranges for different harmonic potentials. The exact values are denoted by the black 
points. The ellipses on the left panel denote the 68% confidence intervals and were determined 
by fitting the ERE to the interpolated phase-shift at a given value of u. The parameters in the 
right panel were determined by fitting the three lowest states of the spectrum to the perturbative 



expressions of sect. II B 



can also be used to extract these parameters. In Fig. [13] the extracted 'scattering volume', 
ai, and 'effective momentum', ri, in the ^Pq channel, determined through the perturbative 
expressions, are shown. The behavior as o;^ — )■ is consistent with the exact result. 
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V. CONCLUSION 



The NN phase-shift below the inelastic threshold can be determined from the eigenvalue 
spectrum of two interacting nucleons confined to move in a harmonic oscillator potential. 
The conventional discussions of scattering from a potential that falls faster than 1/r, and 
the connection between the scattering amplitude and the location of poles in the complex 
energy-plane corresponding to bound-states is complicated by the fact that the harmonic 
potential is confining and asymptotic scattering states cannot be defined for any non-zero 
value of the harmonic oscillator frequency, u. As a result, the zero-range relation between 
the scattering phase-shift and the energy-eigenvalues is modified by the non-zero value of 
the harmonic potential within the range of the nuclear interaction, giving rise to finite- 
range corrections. These corrections are not present when dealing with a pionless EFT 
description of two nucleons confined within a harmonic trap (when the cutoff is taken to 
infinity). However, any pionful theory of the nuclear interaction that describes nuclear 
processes above the t-channel cut will have to address these finite-range issues. 

We have studied these aspects numerically for two nucleons confined by a harmonic po- 
tential. The nuclear interaction was modeled by the JISP16 potential, but our results are 
general and can be applied to other phenomenological or chiral effective field theory inter- 
actions. We have explored uncoupled channels and found that for small values of cj, the 
low-energy phase shift can be extracted from the energy-eigenvalues through an extrapola- 
tion to a; = 0. At the level of precision to which we have performed the calculations, the 
energy-eigenvalues combined with the zero-range relation supplemented by an extrapolation 
to a; = allow for the determination of the low-energy NN elastic scattering phase shifts. 
Further, such calculations enable a precise determination of the deuteron binding energy. 

Since the methods we present here are clearly non-perturbative and include all antisym- 
metrization effects, an interesting application of eq. ^ would be to the elastic scattering 
of two nuclear systems, with one or both composed of more than one nucleon, below in- 
elastic and re-arrangement thresholds [IS]. The processes we have in mind are nd, nt and 
na scattering. Calculations of three-, four-, and five-nucleon systems can be performed 
within harmonic potentials with small u (to access low-energy phase shifts and minimize 
finite-range effects), and an application of eq. (|2]), modified by the reduced mass and an 
appropriate subtraction for the center-of-mass energy, and extrapolation to a; = would 
give the scattering phase-shift at low energies. This method contrasts with those currently 
in use, such as Faddeev [39], Faddeev-Yakubovsky |10j, AGS ^41j, Hyperspherical Harmon- 
ics na Ha], NCSM/RGM m], GFMC [15], and J-matrix methods [16]. There remains 
technical challenges to obtaining sufficient convergence with increasing Nms and/or to ex- 
tending the corrections for finite u to higher order terms. A possible strategy to alleviate 
such issues is through the use of HO-based EFT methods [13], [20]. Work in this direction is 
under way. 
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